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Abstract 

The role of collective phonon excitations on the properties of cold atoms in optical lattices 
is investigated. These phonon excitations are collective excitations, whose appearance is 
caused by intersite atomic interactions correlating the atoms, and they do not arise without 
such interactions. These collective excitations should not be confused with lattice vibrations 
produced by an external force. No such a force is assumed. But the considered phonons are 
purely self-organized collective excitations, characterizing atomic oscillations around lattice 
sites, due to intersite atomic interactions. It is shown that these excitations can essentially 
influence the possibility of atoms to be localized. The states that would be insulating in 
the absence of phonon excitations can become delocalized when these excitations are taken 
into account. This concerns long-range as well as local atomic interactions. To characterize 
the region of stability, the Lindemann criterion is used. 
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1 Introduction 


Cold atoms in optical lattices are usually considered in the frame of the Hubbard model (see, e.g., 
reviews [Ml)- The periodic potential of an optical lattice is imposed by external laser beams. 
This potential is hxed in space, prescribing a lattice formed by the lattice sites a*. We do not 
assume the existence of external helds that would move the lattice. 

If atomic interactions between different sites are neglected, the individual properties of atoms 
in each potential well are completely prescribed by the given optical lattice. However, as soon 
as the intersite atomic interactions are taken into account, the low-energy positions of the atoms 
are not exactly those of the potential minima of the optical lattice. In other words, an atom 
in a potential well experiences an oscillational motion due to the interaction with the other 
atoms. These oscillations, due to the interaction effect, can be characterized as collective phonon 
excitations. 

It is these collective excitations that are considered in our paper. As is well known, collective 
excitations in many cases can essentially influence the system stability. Our aim is to study 
the role of such phonon excitations for atoms in optical lattices. We use the standard method 
of taking into account collective excitations, by considering small deviations from equilibrium 
values. 

It is important to stress the necessity of intersite atomic interactions, without which phonon 
excitations cannot exist. The situations of externally shaking the lattice and of the existence 
of self-organized collective excitations in a system of correlated atoms are principally different 
and should not be confused. In the former case, vibrations should be produced by an external 
held and would exist without any intersite atomic interactions. But in the latter case, there is 
no external shaking held and the collective excitations do not exist in the absence of intersite 
atomic interactions. 

It is the aim of the present paper to study the properties of insulating atomic states in optical 
lattices, taking into account the arising phonon excitations and their inhuence on the stability of 
the insulating states. 

It turns out that phonon excitations can essentially inhuence the properties of atomic states 
in optical lattices. Such excitations play an important role in dehning the boundary of the 
region where atoms can be localized. The presence of phonons can provoke an instability of 
an insulating state, triggering atomic delocalization, often destroying the insulating state that 
would exist without these excitations. This delocalization ehect can occur for both types of 
atomic interactions, long-range as well as local. 

Throughout the paper, the system of units is employed, where the Planck and Boltzmann 
constants equal to one {h = 1, = !)• 

2 Main definitions and notations 

In this section, we give the main dehnitions and notations that will be used in the following 
sections. 

We consider a hxed optical lattice described by the spatial points {aj} corresponding to 
lattice sites enumerated by the index i = 1,2,... ,Nl. The elementary lattice cell is characterized 
by the set of vectors a = {oq}, where the spatial components are enumerated by the index 
a = 1,2,... ,d. For the sake of generality, we consider a d-dimensional space, which will make it 
straightforward to analyze the particular cases of one- two- and three-dimensional lattices. 
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The lattice contains N atoms, whose ratio to the number of sites Nl dehnes the hlling factor 


iy^ — = pa, UJ 

in which p is the average atomic density and a is the mean interatomic distance, given by the 
expressions 


with V being the system volume. The hlling factor can be an arbitrary positive number. 
The optical lattice, formed by laser beams, is characterized by the lattice potential 


VL(r) = ^ Vasin^ {.Kto) , (3) 

a.=l 

where = n/a^. The recoil energy, playing the role of a characteristic kinetic energy, is denoted 


Er = ^ 

^ o_ 


1 “ 

i;2 = lV- 


al 

a=l “ 


with m being atomic mass. The Schrodinger lattice Hamiltonian is 

fflW = -|^ + n(r), (5) 

with the lattice potential (3). 

The extended Hubbard Hamiltonian is derived in the usual way. One starts with the standard 
Hamiltonian of atoms with pair interactions ‘h(r), expands the held operators over the Wannier 
functions w{r — a*), and restricts the consideration to the lowest energy band. For what follows, 
we shall need the notations for the matrix elements over the Wannier functions: the hopping 
term 

Jij = - j - a'i)HL{r)w{T - a^) dr , (6) 

the intersite interaction 


Uij = / |tc(r — aj)f <F(r — r')|t(;(r'— aj)f drdr', 


the momentum squared 


- aj)(-V^)tc(r - a^) dr , 


and the average lattice parameter 


Vl= w*{r)VL{r)w{r) dr 


With these notations, one obtains the extended Hubbard Hamiltonian 


ic]cjCi , 
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in which U = Ujj. The constant term Vl in Hamiltonian (10) can be omitted. In the hopping 
term, it is customary to consider only the nearest neighbors, denoting as J the value of Jtj related 
to these nearest neighbors. Neglecting the last term in Eq. (10), describing intersite atomic 
interactions, would reduce the Hamiltonian to the standard Hubbard model. This omission can 
be motivated by the fact that usually the value of Uij for i j is smaller than the on-site 
interaction U. However, for the treatment of phonon excitations, the intersite interactions are 
crucial, even when they are small. 

The atomic interactions, generally, contain two parts, 

<h(r) = <hzoc(r) + *f>non(r) , (11) 

the local interactions, described by a delta function. 


d);„e(r) = <h,^(r) , (12) 

and non-local long-range interactions <h„o„(r), such as dipolar interactions [SHH]. The parameters 
of the local and non-local interactions can be connected, but to a large extent, the two types of 
the interactions can be treated as independent. Moreover, the strengths of these interactions can 
be varied in a wide range. Respectively, the interaction term (7) is a sum of two parts, 

u., = t/y + c/j”", ( 13 ) 


corresponding to the local interactions, 

Ul°'' = ^d j \w{r - ai)\‘^\w{r - dr, 


and to non-local interactions. 


Tjnon 

^ij 


w{r - ai)p<h„o„(r - r')|tc(r' - drdr' . 


(14) 


(15) 


The strength of the local interaction potential depends on the system dimensionality and 
setup geometry. Considering a system with atomic dynamics in d dimensions, we can keep in 
mind the realistic situation, when 3 — d directions are conhned to the ground state of a harmonic 
oscillator, with a frequency u±, to a size l± = 1/ Then, we will have quasi-one-dimensional 
or quasi-two-dimensional systems iiiio]. 

For example, in three dimensions, the local-potential strength is 

$3 = $0 = dvr — , 
m 

where is the s-wave scattering length. For quasi-two dimensional bosons m, one has 


$0 ^ 




l_L — ttg ln[(27r)T2p/j_a^] 

And for quasi-one-dimensional bosons [12], one gets 

<l>o 


$1 = 


27r/j_(/j_ — O.dbOs 
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If the scattering length is much shorter than the length of the transverse conhnement, then the 
above equations reduce to the formula 


= 


•hr 


tts 


All the cases considered above can be summarized in the form 

^eff 


<^>d = 


in which the effective strength 


«e// 


$e// = dvr 


m 


(16) 


(17) 


is expressed through the effective scattering length. The latter in the quasi-one-dimensional case 
reads as 


«e// - 


1 — 0.46as//j 


(d=l), 


in the quasi-two-dimensional case, it is 

®e// 


1 - (Qs/ l±) ln[(27r)3/2p/_La^ 

and in three dimensions, it reduces to 


(d=2), 


(18) 


(19) 


aeff (d ^ 3) . (20) 

Above, we have listed the expressions that we will need in what follows. The derivation of these 
formulas can be found in the cited literature. Typical expressions for the system parameters for 
an insulating state, such as the hopping term and intersite interactions, are given in the Appendix 

A. 

3 Vibrational collective excitations 

As is well known, collective excitations can essentially influence the system properties, dehning 
the stability boundaries of different physical states PI- Thus, the stability of an insulating state 
can depend on the existence of phonon excitations. These collective excitations can modify the 
properties of many-body systems, even when particle interactions are rather small. In particular, 
phonon excitations can destabilize the system leading to particle delocalization destroying an 
insulating state. As examples of systems where phonon excitations can strongly influence the 
region of stability, we can mention ferroelectrics and atoms in double-well potentials 

[min]. Here, we consider the role of phonons on the stability of insulating states in optical 
lattices. Here and in what follows, speaking about insulating states, we keep in mind the states 
typical of Mott insulators. I such states, atoms are well localized. But the notion of localization is 
more general and essentially depends on atomic interactions. We employ the term ” localization” 
in this general sense. 

Collective excitations are standardly introduced by considering small deviations from equi¬ 
librium values. Taking account of phonon excitations can be done in the traditional way [T81IT9] 
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by considering vibrating atoms characterized by vectors r^, oscillating around the related lattice 
sites a^, with introducing the deviations Uj according to the rule 

r^- = aj + Uj- , (21) 

requiring the validity of the conditions for the averages 


at = (rj) , (uj) = 0 , (22) 

where (...) = Tr[...e~^^]/Tr[e~^^] is an average with respect to the total Hamiltonian including 
all thermal and quantum fluctuations. A vanishing average Uj is based on the condition that 
the system is in stable equilibrium, where is the dehnition of the lattice vectors prescribed 
by the given equilibrium optical lattice. The deviation Uj can become nonzero at the points of 
instability, such as the Peierls instability. 

The quantities Uivij) and J(jij) depend on the difference of the spatial variables 

Yij = Yi- Yj = aij + Uij , (23) 


where we use the notation 

up — iij , a^j a^ aj . 

It is worth mentioning the main difference between the optical lattice, whose periodicity and 
lattice vectors a* are strictly prescribed by the imposed laser beams, and a self-organized crystal, 
whose lattice vectors a* are dehned self-consistently through the minimization of a thermody¬ 
namic potential. 

In a localized state, atomic deviations from the lattice sites are supposed to be small, which 
justihes the expansion of the interaction potentials in powers of the deviations. As usual, re¬ 
stricting such an expansion by the second order, we have 

C/(r.,) c. C/., + 5^ £/««“ - i . 

a aj3 


J {Yij) ~ 4- + ^ 


aj3 a 0 

. .1111.. 


(25) 


«/3 


where we use the notations 
dUr 


f/p = U{a 


Tja _ '^^13 

“ da? 


'ij 

“ da? 




Jij — J 


'tj 


■IJJ ? 


= 


da?da^^ ’ 

^ J 


d^J, 




da?da^ 


As usual, to close the system of equations, it is necessary to decouple the high-order products 
of operators. These higher-order operator products, involving the variables of different nature, 
can be decoupled. Thus, we decouple the atomic and vibrational degrees of freedom involving 
the second-order vibrational variables: 


<jUij4c]cjCi = {u'^jU^j)clc]cjCi + u'^jU^j{clc]cjCi) - {u2ju'^j){clc]cjCi) , 

utjUijclcj = {utjU^dclcj + utjufAclcj) - {ufmfMclcj) , 


p^Acj = {Pj)c]cj + pUclcj) - {pp{c]cj) . 


2/.t. 


J. 


(26) 
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Such a decoupling is motivated by the different physical nature of the atomic and deviation 
operators. Keeping in mind the lattice periodicity, the filling factor can be represented as 


_ N 1 I / t \ 


(27) 


Employing the above decouplings in the Hamiltonian, we meet the combination of terms, for 
which it is convenient to introduce the following notations. Thus, we define the effective hopping 
term 


’Ji] — -Ji] ^ 


exp 


and the effective atomic interactions 


£7«^£7«- 


(28) 


(29) 


ap 


whose valnes are renormalized by the presence of atomic vibrations. These atomic vibrations are 
correlated with each other through the effective interaction matrix 


{c\c]cjCi) - 2Jff{c\cj) 

Atoms produce the effective deformation force 


F" = —r/"rW r- + 2 7“rV- 




(30) 


(31) 


caused by atom-vibration correlations. It is also important to notice that for a function f{sLij), 
depending on the difference SLij, the following properties are valid: 


E 


<9/(a, 


VJ 


i9af 


A 

(9a° 


= 0 , 




da?da^ 


A. 

i9a° 


E 


dfjajj 

(9a£' 


= 0 . 


(32) 


These properties will be used in the final presentation of the system Hamiltonian. 
Accomplishing the described procedure for Hamiltonian (10), we obtain 

H = Fjv + Hat + H-aib + • 

Here the first term is the non-operator quantity 


1 


= i E E *?/ - "Ed); 

i^j a/3 j 

Atoms are described by the renormalized Hamiltonian 


= - E + T E + ■) E UijclcjCjC, + 


(33) 


(34) 


j j 

Collective vibrational degrees of freedom are characterized by the Hamiltonian 


Pi 


(35) 


y.a-t'E2m 4EE*« 


exp ex P 

■ 


(36) 


i^j a(3 
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And the last term Hint corresponds to local deformations caused by the correlations between 
atomic and vibrational degrees of freedom, 

= ( 37 ) 

Hi “ 

Passing from the relative deviations u^- to the single-site deviations u^, and using the above 
properties, the vibrational Hamiltonian part can be represented as 

^ <“7 • (38) 

i Hi q :/3 

The effective deformation force (31) enjoys the property 

F“ = -Fg (* 7^ i) . (39) 

Therefore the deformation term, caused by the correlations between atomic and vibrational 
degrees of freedom, can be rewritten as 

»...< = (40) 

Hi « 

Thus, all terms of the Hamiltonian H are dehned. 

Let us stress that the vibrational collective excitations appear only when there exist intersite 
atomic interactions correlating atoms. In the presence of these interactions, atoms move in an 
effective potential composed of an optical lattice and a self-organized held formed by intersite 
interactions. The optical lattice does not correlate atoms, prescribing only their individual prop¬ 
erties. While intersite atomic interactions do collectivize the atoms, whose collective vibrations 
play the role of collective phonon excitations. 


4 Quantization of phonon variables 


Quantized phonon variables are introduced so that to diagonalize the part of the Hamiltonian 
containing atomic deviations. In our case, the difference with the standard introduction of phonon 
operators is due to the existence of the linear in the deviation term (40). Dealing with such linear 
terms requires to slightly modify the corresponding canonical transformation [20]. In that case, 
the phonon operators are introduced by means of the nonuniform transformation 



in which are the polarization vectors, with s being the polarization index. The phonon 
frequencies are given by the eigenproblem 


u 

m 


HH) 0 


= ^kse 


a 

ks 5 


(42) 
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with the effective interaction matrix (30). 
yields 

A“ = 


Diagonalizing the phonon part of the Hamiltonian 

EE-'t/cf. («) 


ii*i) f 


where we use the notation 

a P 

a/3 _ \ ^ ^ks^ks ik-a^ 

o jy ^ mul 

ks 

and where the effective deformation force, acting on an atom, is 


(44) 


pa 

* o' ■ 


(45) 


The presence of the term Aj in the canonical transformation (41) distinguishes the latter from 
the standard canonical transformation in the quantization of phonon variables. 

It is easy to see that variables (41) satisfy the usual commutation relations 

[<> Pj] = ■ 

Conditions (22) are valid, since 

(A,) = 0 . (46) 

Then Hamiltonian (33) results in the sum 

H = E]\f + Hat + Hph + Hind ■ (47) 

Here the hrst term is the same as in Eq. (34). The atomic Hamiltonian is given by Eq. (35). 
The phonon Hamiltonian is diagonal. 



And the last term is the Hamiltonian of effective multi-atomic interactions, induced by atomic 
vibrations, 

y„.=EE^^S'’d' («) 

aid 

Using properties (32), it is straightforward to check that the average force (45) is zero, that 
is, (F") = 0. Hence, the induced term (49), in the mean-held approximation is zero, {Hind) = 0, 
since 

{FtF^) = {FniFf} = 0 . 

Moreover, the induced term (49) is much smaller compared to the atomic term (35). This implies 
that the induced term (49) does not inhuence much the properties of atoms. 

The properties of phonons depend on their frequency, for which we have the equation 


u 


^ks = 


m 


EE* 


“/3pa Id i’k-a.ij 
ij ^ks^ks^ 


The deviation correlation function becomes 


ks 


^ V- (^ks(^ks 


coth , . 

muks \ 21 / 


(50) 


(51) 
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And the average kinetic energy per atom is 

P'\ 1 

2m j 4z/iV 


^w/i.coth 


(52) 


ks 


An important quantity, characterizing the width of atomic vibrations, is the mean-square 
deviation tq dehned by the equation 

d 


'o — 


= ■ 


(53) 


Q.= l 


This quantity should not be confused with /q that dehnes the width of a wave packet at a lattice 
site. While tq is the mean deviation of this wave packet oscillating around the lattice site. 

In the case of a d-dimensional cubic lattice, the frequency is the same for all polarizations, 
which can be described by the relation 

d 


= 


d 


E 

S=1 


2 

^ks 


The latter results in the equation for the phonon frequency 


Ui. = 


m 


gik-aij 


(54) 


(55) 


with the dynamical matrix 


Dij — 


d 


E* 


OlOL 


Q:=l 


d “ dafdaf 

a=l ^ ^ 


For a cubic lattice, the mean-square deviation (53) reads as 


ud 


^0 = 




dk 


(56) 


(57) 


2mp Je Wfc \2TJ {27Ty ’ 

with the integration over the Brillouin zone. Taking into account only the nearest-neighbor 
interactions leads to the effective phonon dispersion 

k„a' 


9 dZZ V-\ 

^k= - Dq ^ 

OL 


Sin 


(58) 


where Oq, = a and Dq is for the nearest neighbors. In the long-wave limit, Eq. (58) reduces 
to the acoustic spectrum 


Uk ^ Cq/c 




^ 0 


(59) 


with the sound velocity 



(60) 


The value of the mean-square deviation (57) describes the properties of the localized state 
and dehnes the region, where it can exist. As is clear, to be treated as localized, the state has to 
enjoy a mean-square deviation that would be much smaller than the mean interatomic distance. 
This is the essence of the Lindemann criterion of stability that will be considered in the next 
section. 
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5 Possibility of phonon instability 


One of the most important characteristics of a localized solid-like state is the mean-square de¬ 
viation (53) or (57). A localized state can exist only when this deviation rg is much smaller 
than the distance a between the nearest neighbors. This statement is the well known Lindemann 
criterion of lattice stability (see, e.g., [iHiiiniEi]). According to this criterion, the majority of 
solids become unstable and melt when the Lindemann ratio ro/a surpasses 0.2. This criterion is 
valid for anharmonic crystals as well |22] . Even for such a strongly anharmonic quantum crystal 
as ^He the Lindemann ratio is 0.3, which is measured experimentally [23]. In the weakest form, 
the Lindemann criterion states that for the stability of a localized solid-like system it is 
necessary that 

— < 1 . (61) 
a 

The meaning of the Lindemann criterion is evident: If the mean-square deviations of neighbor¬ 
ing atoms would be comparable to their mean interatomic distance, the system could not be 
considered as localized. 

In the long-wave limit, the phonon spectrum is acoustic, as is shown in Eq. (59). This tells 
us that, calculating the mean-square deviation (57), the limit of small wave vectors can produce 
infrared divergence, depending on the system dimensionality and temperature. In order to study 
when and how this happens, we can employ the standard procedure of limiting integral (57) from 
below by introducing the minimal wave vector kmin that is assumed to tend to zero. Equivalently, 
it is possible to dehne the minimal wave vector as kmin = tt/L, where L = is the length of 

the lattice. As is clear, for a large lattice, with the number of lattice sites Nl oo, the minimal 
wave vector tends to zero. 

Considering integral (57) at hnite temperatures T > 0 for low dimensionality d < 2 shows 
that the mean-square deviation diverges as 


'o — 


2^7r2(2-d)z/Do 


(d < 2, T > 0) , 


where Do is the dynamical matrix, as in Eq. (58), and Nl —)■ cxo. In particular, for one¬ 
dimensional space (d = 1), the divergence is linear in the number of sites Nl, 


'o — 


TNl 

2'k'^i'Do 


(d= 1, T > 0) . 


And for two-dimensional space (d = 2), the divergence is logarithmic in Nl, 


(62) 



TXyiNl 

(27r)2z/T)o 


(d = 2, T > 0) . 


(63) 


This means that at hnite temperatures the localized state is unstable in dimensions d = 1 and 
d = 2 for asymptotically large lattices, where Nl —>■ C) 0 . Although the lattice could exist for such 
Nl that would be large, but hnite, at the same time satisfying the Lindemann criterion (61). 

Let us notice that the dynamical matrix Dq, of course, depends on the parameters of the 
optical lattice. This will be shown below by explicit equations. However the infrared divergence, 
considered above, is the feature typical of low-dimensional systems. Recall that, according to the 
Mermin-Wagner theorem [25ll^ continuous symmetry at hnite temperature cannot be broken 
in the spaces of dimensionality lower than three (d < 2), irrespectively of the strength of their 
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interactions, provided that the latter are of short-range type. Such a low-dimensional infrared 
instability is a purely dimensional effect. However, this low-dimensional effect is not of great 
danger for real world that is three-dimensional. Actually, dealing with optical lattices, one 
always deals with three-dimensional systems, which can be reduced to quasi-one-dimensional or 
quasi-two-dimensional by integrating out some degrees of freedom. 

It is generally accepted that for solid states the Debye approximation gives a quite accurate 
description [T8 | IT ^ I2T1I23] . In this approximation, the integration over the Brillouin zone is 
replaced by the integration over the Debye sphere. 



2 

{4n)‘‘/^r{d/2) 


(64) 


limited by the Debye radius ko that is dehned by the normalization condition 

f dk Nl P 


V V ’ 


where p is the average density 

This gives the Debye radius 


N V 
V ^ 


P=T7= . 


ko — 


\/4^ 


irii 


l/d 


The spectrum is taken to be isotropic, with the frequency 

cjfc = Cok (0 < fc < ko) , 
whose upper limit dehnes the Debye temperature 

Td = Coko ■ 


(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 


With the notation for the Debye radius (67), the replacement (64) takes the form 

r dk d 


(27r) 


-)■ 


k'^-^ dk . 


(koa)^ 


(70) 


Below, we study in more detail the Lindemann criterion, employing the Debye approximation. 
At zero temperature, the mean-square deviation becomes 




^0 = 


2{d — l)mTo 
The Lindemann criterion of stability (61) yields 

d^ 


(T = 0) . 


(71) 


T„> 


2{d — l)ma? 


(r = o) 


(72) 


The meaning of the latter inequality is very transparent: a localized state can be formed only 
when the effective potential energy is larger than the kinetic energy of atoms. As Eq. (72) 
demonstrates, no localized state can exist at zero temperature for <7=1. 
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At high temperatures, the mean-square deviation is given by the expression 


I n — 


T(f 


° {d - 2)mTl 
Then the Lindemann stability criterion gives 


Tn > 


Td2 


{d — 2)ma? 


(T > Td) . 


( 73 ) 


{T » Td) . 


(74) 


This tells us that, at such temperatures, there can be no localized state for d = 2. At these 
temperatures, only a three-dimensional localized state can exist, provided that 


Tn > 


9T 


ma^ 


(T > Td, d — 3) 


(75) 


To simplify the consideration, for a well localized insulating state, one can use the averages 


(clcj) = SijU and 


{clcUjd) = {clci){c]cj) = u^ j) . 


As a result, the interaction matrix (30) reduces to 


(76) 


(77) 


An interesting question is how atomic vibrations influence a localized state due to local 
interactions. For such a case, considering a cubic lattice, taking account of only the nearest 
neighbors, and involving the formulas from the Appendix A, we get 


Uij = U exp ( — 


a^d 

af 


(78) 


where the equality = a^d is used and U = Uioc is defined in the Appendix A. The dynamical 
matrix (56) becomes 


ua\ ( a^d 

^0 = ^ t/ exp - ^ 


The sound velocity (60) reads as 


Co = 


va" 


V 

m 


211 


o?d 


I - ill 


(79) 


(80) 


which can also be represented in the form 


4z/Ja2 


Co = 


(tt^ — 4)Vo/od V ^ 


-U, 


( 81 ) 


through the parameters of the optical lattice. 

For a cubic lattice at zero temperature, the Lindemann criterion of stability (61) can be 
written as 


16(27r)^/^z/^/^ J 
(7r2-4)K)ci 


d 


1/1 

h 


®e// 


j3—d id—2 

i_|_ tg 


> 


d^ 


d-l 


(82) 
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Consequently, the system stability essentially depends on the lattice parameters, space dimen¬ 
sionality, as well as on atomic interactions. 

At zero temperature, there can arise only two- and three-dimensional localized states. Con¬ 
sidering the two-dimensional case {d = 2), we have 


$2 — 

U = 


^eff 

l± ’ 
^eff 

(27r)3/2/^/2 


The Debye temperature (69) becomes 


\/4^ 

kn — - 

a 


(d = 2). 


... Si^’Ar 

The Lindemann criterion (72) is equivalent to the inequality 

ma^To > 2 , 



(83) 


(84) 


which gives 

Comparing the ratio Uioc/J, given in the Appendix A, with criterion (85), we see that the 
atomic state can be insulating, when no phonon degrees of freedom are taken into account and 
U is much larger than J. But as soon as phonon excitations are included, the stability condition 
(85), for the same parameters, may be not valid, which means that atoms delocalize. 

For the three-dimensional case (d = 3), we get 


U = 


The Debye temperature (69) becomes 


*^3 — *^o ; kn — 
$0 


(27r)3/2/3 


(67r^)^/3 


(d = 3) . 


va va 


Tn = 3.482 ^ / — exp l - 


3a^ 

HI 


mlg V Iq 

At zero temperature, the Lindemann criterion (72) takes the form 


( 86 ) 


(87) 


ma^Ti) > 


which yields 


u 


exp > 0.648. 


9 

4 ’ 


3a^ 


Hi 


I I 4 72 

‘o/ V D 

Again, the system can be insulating without phonons, but becoming delocalized in the presence 
of the latter. 
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Using the formulas of the Appendix A, the localization conditions (85) and (88) can be 
represented in another form by taking into account the expressions for the hopping rate 

J = 2.935 Vo exp {d = 2) , 

J = 4.402 Uo exp (d = 3) . 

Then Eq. (85) yields 

while criterion (88) gives 

a 8.137 /Uo\Y/oV 

This shows that the Lindemann criterion of stability requires that atomic interactions, hence the 
scattering length, be sufficiently large. If these conditions are not valid, the insulating state can 
be destroyed by phonon vibrations. 

In order to better understand the physics of the phonon instability, we need to analyze how the 
occurrence of atomic vibrations influences the values of the hopping parameter and interaction 
matrix. 




6 Renormalization of atomic parameters 


Phonon excitations renormalize atomic parameters according to Eqs. (28) and (29). The renor¬ 
malized quantities are shifted, because of the phonon existence, resulting in 

Jij = Jij + AJjj , Uij = Uij + ^Uij , (91) 


where, taking into account that 

«Y) = 2(1 - , 

we have 

= - E -^7 . At/y = - . 

0/3 0/3 

In the Debye approximation, we hnd 

which gives 

2 2 

AT-—- — AU- • = - ^ /■/"" 

“ d ^ d ^ ' 


In the tight-binding approximation, where /q a, the hopping parameter shift is 

^ 2^2 

4/4 ’ 


(92) 


(93) 

(94) 


(95) 
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which shows that the hopping parameter increases due to collective phonon excitations. At zero 
temperature, the relative shift is 


AJ, 






Jij 8{d — l)m/oTD 

For a two-dimensional lattice {d = 2), this gives 


(96) 


A Jij „ a 

^ = 0.079 


/ 


J, 




z//o V uaeff V2/0 


= 0.074 


aVo 


Z/3/2/oJ Y aeff ’ 


(97) 


and for a three-dimensional lattice (d = 3), the relative shift is 


^4' nnr:. ® MO / 30^ ^ aid) jk 

—— = 0.054 — W- exp ^ = 0.237 , \ — 

Jij ulo V V4io/ V Os 


(98) 


The increase of the hopping parameter can be rather noticeable. Thus, for tq/Iq ~ 0.3 and 
/o/a ~ 0.1, shift (95) is of the order of Jij. 

The shift of the interaction matrix, with the local interaction potential, reads as 


AU,, 


2.„2 


a r, 



4 




which means that interaction matrix increases. At zero temperature, this yields 


(99) 


AUij aJdJ 

Uij 2{d — l)mlQTD 


( 100 ) 


The increase of the effective interaction matrix Uij, caused by collective phonon excitations, 
can be sufficiently large, even in the case of the local atomic interactions. For instance, if 
ro/k 0.3 and k/a ~ 0.1, then the shift (99) is of the order of Uij. Thus, the phonon vibrations 
can essentially renormalize the atomic Hamiltonian parameters. A simple example of temperature 
dependence is given in Appendix B. 

In this way, collective phonon excitations lead to the increase of the hopping parameter and 
of the interaction matrix corresponding to the interactions of atoms in different lattice sites. 
However, notice that the on-site atomic interaction parameter U in Hamiltonian (35) remains 
unchanged. 


7 Physics of phonon instability 

Now it is straightforward to understand why collective phonon excitations can lead to the in¬ 
stability of an insulating state. To this end, keeping in mind low temperatures, let us compare 
the energy of the system described by the Hamiltonian without phonon degrees of freedom with 
the energy of the system including phonon excitations. These energies are given by the average 
values of the related Hamiltonians. The energy of the system with phonons is dehned as the 
average E = (77) of Hamiltonian (33). While the energy E of the system without phonons can 
be dehned as the average of Hamiltonian (10). So, we need to consider the difference A = E — E. 
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As is explained above, when analyzing the interaction Hamiltonian (40), we find that its 
average value is small, {Hint) ~ 0. The Appendix A shows that the intersite term 17^, with i ^ j 
is much smaller than the on-site term U. Thus, we obtain the energy difference 

AE ^ AJij{clcj) . 

Hi 

In the previous section, it is shown that the hopping parameter shift is positive. Thus, the 
existence of phonons increases the hopping parameter and, consequently, decreases the system 
energy. Since the system with phonons prefers a smaller energy, it is more stable than the system 
without these excitations. 

A localized state is metastable, since in order to destroy it, atoms have to penetrate through 
the barrier Vq created by the optical lattice. Such a penetration requires some time that charac¬ 
terizes the lifetime of the metastable state. This lifetime can be estimated [28H32] as 

, [Vb 

tract = To exp - 

\^N 

where Vg is the barrier height, is the characteristic energy of noise, which is defined by the 
characteristic kinetic energy, and tq is the period of atomic oscillations at the bottom of the 
well. For classical systems the characteristic noise energy Sj^ coincides with temperature T, 
which results in the Arrhenius formula. For quantum systems, the characteristic kinetic energy 
is defined by the energy that, in our case, is the recoil energy Er. And the barrier height for an 
optical lattice is Vq. With the oscillation period r = 27r/a;o, we have the lifetime 

271 

tract = - exp 

UJo 

As an illustration of the above consideration, let us study some experiments, where the 
insulating state in an optical lattice has been observed. For example, let us consider the exper¬ 
iment of Greiner et al. [33], where a three-dimensional cubic lattice was created, loaded with 
the atoms of ®^Rb, with the filling number close to one. Varying the optical-lattice parameters, 
both the localized insulating and delocalized superfiuid states were realized. The lattice was 
formed by laser beams of the wavelength A = 0.852 x 10“^ cm, which makes the lattice param¬ 
eter a = A/2 = 0.426 x 10“^ cm. With the mass m = 1.443 x 10“^^ g, the recoil energy is 
Er = 2.093 X 10“^^ erg. At the lattice depth Vq = ISEr = 2.721 x 10“^^ erg, atoms are local¬ 
ized in a Mott state. Under these parameters, the effective frequency uq = 1.509 x 10“^^ erg, 
which gives the wave packet width Iq = 0.715 x 10“^ cm. Since lo/a = 0.168, the tight-binding 
approximation is applicable. With the scattering length = 0.545 x 10“® cm, we have the ratio 
ttg/lo = 0.076. 

At the same time, for Vq = ISEr numerical calculations [3S| give U = 0.35i7g and J = 
77/36 = 0.97 X 10~^Er. Then the Lindemann criterion (90) requires that ag/lo be larger than 
330. According to this criterion, the localized insulating state is unstable if the phonon degrees of 
freedom are taken into account. Therefore, the experimentally observed localized states of atoms 
interacting through local delta-function forces, are not absolutely stable equilibrium states, but 
rather are metastable states. 

The oscillation time of an atom in a well is r = 27r/a;o = 0.44 x 10“^ s. Therefore, the lifetime 
of a metastable insulating state, under Vq = 13Er is ~ 200 s. This is a quite long time, 
allowing for its easy observation. Moreover, this time is longer than the typical lifetime of atoms 
in a trap, which is of order of seconds or tens of seconds [MIISS] . Therefore, it is very feasible to 
create long-lived metastable insulating states of trapped atoms in optical lattices. 
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8 Atoms with dipole interactions 


In the case of non-local long-range atomic interactions <h„o„(r), such as dipolar interactions, 
atoms in an optical lattice, where Iq <C a, can be described [SHH] by the extended Hubbard model 
with the effective interactions 

U^J = % ( 101 ) 




in which Cd oc /Iq, with /tq being the magnetic (or electric) atomic dipole. Since dipole interac¬ 
tions are of hard-core type, each lattice site can host only a single atom, that is, the hlling factor 
is one, p = 1. For the interactions of form (101), one has 


t/“ = -?,Cd 


5 ’ 

ij 


= 


SCd ( 


O \ V 

The shift of the effective interactions, due to phonon excitations, according to Eq. (94), reads as 


3(6-d) ,Co 
^ <i “ al 


( 102 ) 


Taking into account only the nearest neighbors, yields 


AHi,- = 




(103) 


Again, we see that the presence of phonons increases effective interactions. This increase is not 
large, since, for vq a, Eq. (103) shows that the shift is much smaller than Uij. 

For a well localized insulating state, the interaction matrix (30) leads to 








(104) 


which defines the dynamical matrix (56) as 


Dij — 


3(5 — d) 


d 




Taking into account only the nearest neighbors gives 

3(5-d) Cd 


T)n = 


d7/2 


The sound velocity (60) becomes 


Co = 


3(5 - d)CD 

ma^d’^C 


1/2 


which for the Debye temperature (69) yields 


Td = 


2^37r(5 - d) 

dJC 


d^(d 

2 I 2 


l/d 




D 


ma° 


(105) 


(106) 


(107) 
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In particular, for two and three dimensions, we find 


Tz) = 3.162 {d = 2), 

V ma^ 

Tn = 1.396 id = 3). 

V ma^ 

It is convenient to introduce the dipole length 

mCn 


( 108 ) 


(109) 


which sometimes is also called the effective dipole scattering length. Then the Lindemann cri¬ 
terion (61) for the stability of a localized state can be written, depending on dimensionality, 

(<i = 2), 

(d = i). ( 110 ) 

This tells us that the effective atomic interaction, i.e., the dipole length, has to be sufficiently 
strong for the stability of the localized state. 

In order to estimate the characteristic interaction parameters, let us consider the systems of 
cold trapped atoms of ^^Cr [36], ^®^Er [37], and ^®^Dy [38]. We keep in mind that the typical 
nearest-neighbor distance in an optical lattice is a 10 ® cm. 

The dipole magnetic moment of ^^Cr is po = 6 / 15 , which gives po = 5.564 x 10“^° erg/G. 
The atomic mass is m = 0.863 x 10“^^ g. The scattering length is = 103 as- The sound 
velocity is cq ~ 0.2 cm/s. The Debye frequency is ud ~ 0.6 x 10® 1/s. The Debye temperature 
is To ~ 0.7 X 10“^^ erg or To ~ 0.5 x 10“® K. Taking Co ~ we get ao = 2.4 x 10“'^ cm. This 
is much shorter than the typical lattice distance a. That is, the stability criterion (110) cannot 
be satished. 

For ^®®Er, the dipole moment is po = 7ps, which yields po = 6.492 x 10“^*^ erg/G. The mass 
is m = 2.777 x 10“^^ g. Then the dipole length is a^) = 1.052 x 10“® cm. This is also shorter 
than the typical intersite distance a. Hence the stability criterion (110) is again not valid. 

In the case of ^^'^Dy, the dipole moment is po = 10/iB, which results in po = 9.274 x 10“^° 
erg/G. The mass is m = 2.698 x 10“^^ g. This gives the dipole length ao = 2.087 x 10“® cm that 
again does not satisfy the stability criterion ( 110 ). 

Thus, phonon excitations do not allow for the formation of localized states for the above 
atoms with pure dipolar forces. But there exist polar molecules for which magnetic (or electric) 
moments can reach lOO/is [7]. Such polar molecules, with the dipolar lengths several orders 
larger than those of the above atoms, can satisfy the stability criterion, hence, can form localized 
states that are stable against phonon collective excitations. And even when absolutely stable 
states are not allowed, there can exist very long-lived metastable states. 


^>0.4 

a 

ao 


> 2.6 


a 


9 Conclusion 

We have studied the influence of phonon collective excitations on the possibility of cold atoms 
to form localized states in optical lattices. The appearing phonon excitations are self-organized 
atomic fluctuations, caused by intersite atomic interactions collectivizing the atoms. 
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It turns out that such phonon excitations are very important even for atoms with local 
delta-function interactions. Phonon oscillations can destroy the insulating state that would exist 
without them. 

The physical mechanism, by which phonons can destabilize an insulating state, is the fact 
that phonons decrease the system energy by increasing the hopping parameter. 

The localized state can be stabilized by strong atomic interactions. These conclusions are valid 
for both local as well as long-range dipolar interactions. Taking into account collective phonon 
excitations is necessary when studying whether cold atoms in optical lattices can form localized 
insulating states in equilibrium. Even when in absolute equilibrium a localized state cannot 
exist, being unstable with respect to phonon excitations, this does not prohibit the existence of 
metastable localized states in optical lattices, which live so long that they can be easily observed 
and studied. 

The reason why the conditions for the existence of localized states for a pure Hubbard model 
without phonons and that one with the phonon excitations are different is easy to understand. 
In a pure Hubbard model, the state is localized when the width Iq of the wave packet at a lattice 
site is narrow, such that it is much smaller than the intersite distance a, or when fco^o ^ 1- The 
latter, since fco^o = (Er/Vq)^^'^, implies that Vq 3> Er. Exactly this condition is observed in the 
experiments studying the insulating Mott state in optical lattices [T1H1I33]. 

However, the wave packet can be narrow, but strongly oscillating around a lattice site. Such 
oscillations are characterized by the mean-square deviation tq. Then the localized state can exist 
if these oscillating wave packets do not intersect with each other, which is the meaning of the 
Lindemann criterion of stability, requiring that rg/a be at least smaller than 1. The latter imposes 
the condition that atomic interactions be sufficiently strong, satisfying inequalities (89) or (90), 
depending on the space dimensionality. The condition lo a does not directly involve the 
parameters of intersite atomic interactions, while the Lindemann criterion of instability strongly 
depends on the parameters of such interactions. This is why these criteria are principally different 
and do not need to coincide. 
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Appendix A. Hopping term and intersite interactions 

In deriving the extended Hubbard Hamiltonian, one employs the basis of Wannier functions. 
It is known that this basis is convenient because Wannier functions can be made well localized [391 - 
E], such that the tight-binding approximation becomes applicable, when the Wannier functions 
are close to harmonic wave packets. 


tc r = 


]^M;„(r), Wa{r) = exp 


a=l 


m 9 


The wave packets are well localized in the sense that 


— < 1 


Irv — 


y/muTa 


or, in other words, that 

k^la C 1 (^0 = Tt) . 

The value la plays the role of the localization length, or width, for the wave packet. The 
corresponding oscillator frequency Ua can be found by minimizing the energy of the whole atomic 
system |3] , which takes into account atomic interactions. When the packet width is mainly dehned 
by the optical potential, then 


UJa = \ — Va /Cq 
m 


= — 


In this approximation, the hopping term (6) is 




Q;=l 


^(^1 -K 


exp 


lE ^ 


a=l 


and the matrix element (14) becomes 


TjloC _ 

“ (2^ 


n 

ya=l 


— exp 


^E{'¥ 


Q;=l 


with the notation 


r Ck 4 2 I 12 '' / Qi \ 2 

= aj — Uj = {cijj} , cijj = |ajj| = / , 


Q=1 


For a cubic lattice, with Va = Vq, Ua = ujq, and la = lo, the hopping term simplihes to 


^ 8/2 


VqcL ) exp 


4/2 ' ’ 


while the above intersite interaction becomes 


ulr = f/^occxp (- 
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with the on-site interaction, due to the local potential, being 


Uloc 


{2tiY/X ' 


For a cubic lattice, with the frequency Wq dehned by the optical potential, we have 

uJo = 2^/\^=-J-Vo 
a\m 


with Er = 7r^/2ma^. In the case of the nearest neighbors, using the relations 


811 4 ttV^o ’ II S Er 


we hud the hopping term 


and the on-site interaction, caused by the local potential. 


J = Vod{—-1] exp 


Ulnr 


JP ^c^/4 _ ‘^^OO-eff 


The ratio of the on-site interaction to the hopping term reads as 


Uloc _ 4<l>rfexp(a^d/4/g) 

J (tt^ — 4)(27r)^/^ /q Vod ’ 


or it may be presented as 


J - A)d \2a‘^Vy 


exp I — v^mPo 


This shows that, for a sufficiently deep lattice and large scattering length, the value of Uioc can 
be made much larger than J, so that the system would be in a well localized insulating state. 
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Appendix B. Mean-field illustration of phonon instability 


The conventional approach to determine the properties of both, the phonons and atoms, is 
based on a self-consistent evalnation of the self energy (Migdal approximation) [12]. The latter 
provides an effective (or renormalized) energy and its imaginary part, an effective scattering rate. 
Such a static approximation might be insufficient, since it does not take into account thermal 
ffuctuations. It is possible to treat quantum and thermal ffuctuations separately. To this end, we 
can replace the phonon operators bks, &L ^>7 their quantum average: b^s ~ {bks) = v^s and 

~ (^L) = '^ks- approximation, we can keep thermal ffuctuations but ignore quantum 

ffuctuations of the phonons. The atoms, on the other hand, are studied in full quantum dynamics. 
This reduces the grand-canonical ensemble at inverse temperature [3 = l/Zc^T, defined by the 
generating function Tre~^^, to a functional integral with respect to thermal ffuctuations of the 
lattice distortions Uj and a trace with respect to the quantum states of the atoms [13] . 

In order to illustrate how the phonon instability can develop, we consider a simplified model to 
show how phonon ffuctuations rise with temperature. For this purpose we choose as the effective 
quasi-atom Hamiltonian h a hopping term and a term that describes the displacement of atoms 
on nearest-neighbor sites. Assuming a bipartite lattice, the atomic Hamiltonian Hat of Efi- (35) 
can be reduced to the effective form 

h = L + HL)/2, h., = i(F„, - Hl)l2 , 


where is an effective chemical potential of atoms. Within the mean-field approximation, this 
leads to a spatially uniform displacement field u entering the action 


^ + 


In 


1 + - 2e-^^ cosh U^Jhl + hl 


dk 

{2'kY ’ 


with the phonon Hamiltonian of Eq. (38). The calculation of the phonon free energy —/31nZ 
requires the integration over the atomic displacement field u, which can be performed in saddle- 
point approximation fixing the displacement field by the saddle-point equation 


dS _ dHyib 
du°' 9 m " 


sinh -|- ^ 2 ^ hidhi/du^ + h 2 dh 2 /du°‘ dk 

cosh(/9/i) — cosh {j3\/h\ + \/h\ + h\ (27r)'^ 


Phonon ffuctuations around the saddle point are described by the fluctuation matrix 


du°‘du^ 


cosh ((3\/hl + cosh(/3/i) 


cosh(/9/x) — cosh + 


1 2 


{hidhi/du°‘ + h2dh2/du'^){hidhi/du^ -f h2dh2/du^) dk 
^ hl + hl (27r)'^ ■ 

The effect of the phonon fluctuations increases with decreasing eigenvalues of the fluctuation 
matrix. The negative sign in front of the positive integral reflects the fact that the phonon 
fluctuations increase with increasing temperature and can eventually lead to an instability of the 
atomic system. Thus, the fluctuation matrix is related to the Lindemann criterion of Sect. |5l 
This behavior is depicted by the temperature dependence of the integrand in Fig. 1. Eigenvalues 
can even become negative, which indicates a phase transition. The latter could either be a 
structural phase transition of the lattice [^IT5] or a melting transition [2l] . 
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Figure Caption 


The integrand of the last equation of Appendix B as a function of temperature T (in units of 
the hopping energy), describing the increase of phonon fluctuations with rising temperature, for 
+ hi = 0.01, 0.8,1 (from bottom to top), d = 3, and fi = 1.2. 
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Figure 1: The integrand of the last equation of Appendix B as a function of temperature T (in 
units of the hopping energy), describing the increase of phonon fluctuations with rising temper¬ 
ature, for \/h\ + h 2 = 0.01, 0.8,1 (from bottom to top), d = 3, and p = 1.2. 
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